Effect of the reservoir size on gas adsorption in inhomogeneous porous media 
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We study the influence of the relative size of the reservoir on the adsorption isotherms of a fluid 
in disordered or inhomogeneous mesoporous solids. We consider both an atomistic model of a fluid 
in a simple, yet structured pore, whose adsorption isotherms are computed by molecular simulation, 
and a coarse-grained model for adsorption in a disordered mesoporous material, studied by a density 
functional approach in a local mean-field approximation. In both cases, the fiuid inside the porous 
solid exchanges matter with a reservoir of gas that is at the same temperature and chemical potential 
and whose relative size can be varied, and the control parameter is the total number of molecules 
present in the porous sample and in the reservoir. Varying the relative sizes of the reservoir and 
the sample may change the shape of the hysteretic isotherms, leading to a "reentrant" behavior 
compared to the grand-canonical isotherm when the latter displays a jump in density. We relate 
these phenomena to the organization of the metastable states that are accessible for the adsorbed 
fluid at a given chemical potential or density. 



I. INTRODUCTION 

Under confinement in disordered mesoporous materi- 
als, the characteristic time scale for relaxation of a fluid 
can become extremely long. As a result, and although 
not always appreciated, equilibrium is often not attained 
and, accordingly, bona fide thermodynamic transitions 
such as the liquid-gas transition in one-component fluids 
or macroscopic phase separation in mixtures, are unob- 
servable. "Capillary condensation" in disordered solids is 
an out-of-equilibrium phenomenon, as illustrated by the 
irreversibility and hysteresis effects found in experiments. 
One typically observes a hysteresis loop that describes 
the isothermal evolution of the amount of fluid adsorbed 
in the porous solid as a function of the applied pressure, 
with branches that differ on adsorption (filling) and on 
dcsorption (draining). This hysteresis loop appears rate- 
independent, and its size and shape vary with tempera- 
ture as well as with the characteristics of the solid {e.g. 
its porosity) or those of the solid-fluid interaction poten- 
tial. Such a phenomenon is related to the existence of a 
large number of "metastable states" in which the system 
can be trapped on the experimental time scale; evolution 
from one metastable state to another then only occurs as 
a result of the action of the applied pressure (or equiva- 
lently, chemical potential) and it proceeds through a se- 
quence of irreversible cooperative condensation (or evap- 
oration) events, generically denoted as "avalanches" . [ij 
The fact that the location and the shape of the hysteresis 
loop are reproducible in experiments indicates that the 
observation time is smaller than the time to reach the 
global equilibrium state, but is larger than local equi- 
libration processes by which the system settles in one 
metastable state. As a result, the behavior of a fluid 
during filling or draining in disordered mesoporous ma- 
terials can be rationalized by envisaging the evolution of 
the system in a free-energy landscape characterized by 



many local minima, i.e. metastable states. 0, S Q 

The above picture of gas adsorption in disordered 
porous media brings in a strong analogy with the out- 
of-equilibrium response of systems driven by an exter- 
nal force in the presence of impurities or other types 
of quenched disorder. This is for example the case of 
magnetization cycles in ferromagnetic materials when a 
magnetic field is ramped up and down and of hysteretic 
martensitic transformations in alloys; Q in both exam- 
ples, avalanches can be detected through some "crack- 
ling noise", magnetic Barkhausen noise in the former, 
acoustic emission in the latter. In such driven disordered 
systems, one expects the occurrence of out-of-equilibrium 
phase transitions as one changes, on top of the driving 
force, some external parameters such as the temperature 
or the characteristics of the intrinsic disorder {e.g. the 
porosity in a porous solid). Q The branches of the hys- 
teresis loop, in particular the adsorption and the dcsorp- 
tion isotherms, may then display jumps (discontinuities): 
indications for such behavior are for instance seen in the 
adsorption of helium in very light aerogels, [sl lol. [Tol| 

Such discontinuities and out-of-equilibrium phase tran- 
sitions, however, arc theoretically predicted on the ba- 
sis of a grand-canonical set-up (for gas adsorption) in 
which the gas reservoir is infinite. In real experiments, 
the reservoir has a finite size which may not always be 
large enough for considering that the fluid inside the 
porous material is in a grand-canonical ensemble with 
fixed chemical potential. What should be expected in 
such situations ? The evolution of the system among 
metastable states may depend on the specific experimen- 
tal set-up. As far as we know, there has been no sys- 
tematic experimental study of the effect of changing the 
size of the gas reservoir (relative to that of the porous 
medium) and no estimate of the condition under which a 
grand-canonical situation is approximately reached. An- 
swering these questions is the primary goal of the present 
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paper. 

The question of the dependence of the hysteretic re- 
sponse on the chosen control parameter has been a httle 
more studied in other driven disordered systems. Two 
extreme cases have been considered: the response of an 
extensive quantity to a change in the (conjugate) "force" 
and the response of an intensive quantity (a force) to a 
change in the conjugate extensive quantity; the former is 
akin to the grand-canonical situation for a fluid in which 
the chemical potential (or actually, the pressure in the 
gas reservoir) is controlled and the latter to a canonical 
situation in which the number of adsorbed fluid molecules 
is controlled. Examples of such studies are found in the 
context of martensitic transformations in which either 



stress or strain is controlled and in that of Barkhausen 
noise in which either magnetic field or, via some feed-back 
mechanism, magnetization is controlled. Q In all cases, 
the loop obtained with the extensive variable as control 
parameter appears as reentrant when compared to that 
obtained with the force as control parameter. (For a the- 
oretical study, see Ref. [l3|) 

One anticipates that the influence of the relative size 
of the reservoir on the adsorption isotherms of a 
fluid in a disordered porous material is intimately con- 
nected to the organization of the metastable states in 
the adsorbed-density/chemical-potential plane. Indeed, 
when the isotherms are smooth both on adsorption and 
on desorption, one may experimentally probe metastable 
states located inside the main hysteresis loop by study- 
ing ascending and descending "scanning curves" , which 
involve partial filling or draining. Q These curves lead to 
a variety of hysteretic subloops that provide direct ev- 
idence for the presence of metastable states inside the 
main loop. [3, iTsl . Actually, metastable states are 
expected everywhere inside the latter. The situation is 
quite different, however, when the adsorption or the des- 
orption branch displays a jump (in a grand-canonical set- 
ting with a very large reservoir) . A discontinuity prevents 
the realization of scanning curves in some portion of the 
isotherm. As a result, part of the interior of the main 
hysteresis loop is now inaccessible. 

In the present work we address the above questions, 
namely the effect of the relative size of the gas reser- 
voir on the adsorption isotherms of a fluid in a dis- 
ordered or inhomogeous porous solid and the connec- 
tion to the distribution of metastable states inside the 
hysteresis loop. We show that even when the grand- 
canonical isotherms display discontinuous jumps, there 
are branches of metastable states inside the main hys- 
teresis loop and that many of these states can be reached 
by varying the relative size of the reservoir, the smaller 
the reservoir the larger the extent to which the branches 
of metastable states are probed. In particular, the ex- 
tent is maximal when the reservoir is so small that the 
fluid inside the porous solid behaves as in a canonical 
situation of fixed number of adsorbed molecules. When 
compared to the loop obtained for an infinite reservoir 
(grand-canonical situation) and plotted as the amount 



adsorbed versus the chemical potential, the hysteresis 
loop for a finite reservoir size then appears as reentrant. 
On the contrary, when the isotherms are smooth (contin- 
uous), there is no infiucnce of the size of the reservoir. 

The rest of the paper is organized as follows: 

In section II, we introduce the two models that have 
been investigated and give some details on the methods 
used to compute the adsorption isotherms. Wc have con- 
sidered: (i) an atomistic model of a fluid in simple, yet 
structured pore, whose adsorption isotherms are com- 
puted by molecular simulation, and (ii) a coarse-grained 
model for adsorption in a disordered mesoporous mate- 
rial, studied by a density functional approach in a local 
mean-field approximation. In both cases, the fluid inside 
the porous solid exchanges matter with a reservoir of gas 
that is at the same temperature and chemical potential 
and whose relative size can be varied. The overall system 
composed of the sample plus the reservoir is taken in the 
canonical ensemble, with the total number of molecules 
as control parameter. 

In section III, we present the results for the atomistic 
model. The simplicity of the system allows one to get a 
clear interpretation of the physical nature of the branches 
of metastable states, whose number is limited. The exis- 
tence of metastable states is a direct consequence of the 
intrinsic inhomogeneity of the pore space. As the relative 
size of the gas reservoir is varied, we find that the explo- 
ration of the branches changes, the isotherm displaying 
larger jumps in the adsorbed density as the size increases. 

Section IV is devoted to the coarse-grained model of a 
disordered porous material. The number of metastable 
states is now enormous, since it is likely to increase expo- 
nentially with the size of the sample. We find a drastic 
change of behavior between the "strong-disorder" regime 
for which the isotherm is continuous and the "weak- 
disorder" one for which the isotherm, here the adsorption 
isotherm on which we focus our study, displays a jump 
corresponding to a macroscopic avalanche. In the former 
regime, the isotherm docs not display any dependence 
on the relative size of the reservoir whereas in the latter, 
it is reentrant for finite reservoir sizes, the more so as 
one decreases the reservoir size, and it shows increasing 
jumps in adsorbed density with increasing reservoir size. 

Finally, wc summarize our main results and give our 
conclusions in section V. In particular, we discuss the 
relevance of our study to experimental situations and we 
stress the important role of the intrinsic inhomogeneity 
induced by the solid matrix. 



II. MODELS AND METHODS 

A. the set-up: sample plus reservoir 

The principles of the method are the following. One 
considers a starting situation where N particles are 
placed inside two cells which are in thermal equilibrium 
with an infinite heat bath at temperature T. One of 
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the cells represents a porous solid sample of volume Vp 
and the other a reservoir of volume Vr. The total vol- 
ume is then equal to V = Vr + Vp. We allow mass 
exchange between the cells whose volumes are kept fixed 
so that at (exchange) equilibrium the chemical poten- 
tial fip of particles inside the reservoir is equal to the 
chemical potential fj.p of particles inside the porous solid: 
fi = ^p = [ip. We measure the average amount of par- 
ticles Np present in the porous solid. Varying the to- 
tal number of particles TV in small steps, one changes 
both the chemical potential of the cells, /i(iV), and the 
number of adsorbed particles, Np{N). We then moni- 
tor the adsorption isotherm p = Np/Vp (for simplicity, 
we omit the subscript P in the fluid density adsorbed 
in the porous sample) as a function of ^. Introducing 
the ratio a = Vp/V and the average reservoir density 
Pp = Np/Vp, where Np ^ N — Np, one has to satisfy 
for each p the constraint pp = (N/V) ^ ap + (I — a)pp 
with the condition p = pp{pp,T) — pp{p,T). Taking 
the limit a — > corresponds to the grand-canonical en- 
semble for the adsorbed fluid (with a controlled chemical 
potential) whereas taking the limit to 1 corresponds to 
the canonical ensemble for the adsorbed fluid (with a 
controlled number of adsorbed particles). Between these 
two extremes, one has a mixed ensemble, with a reservoir 
of variable size compared to the sample, somewhat remi- 
niscent of the mesoscopic canonical ensemble considered 
by Neimark and coworkers. [l7j More details about the 
algorithms used in our two models are given below. 



B. Atomistic model of a fluid in a single structured 

pore 

An atomistic model of confined fluid was considered 
because it allows a realistic description of the various 
(metastable) states the system can adopt and, possi- 
bly, of the transition mechanisms at a molecular scale. 
The chosen system is a simple atomic fluid (Lennard- 
Jones like) confined in a mesoporous substrate. (The 
model more specifically corresponds to Arg on adsorbed 
in nanoporous solid carbon dioxide. ) p^. [l9l| A fully real- 
istic model of the substrate should take into account the 
surface roughness, pore morphology (pore size distribu- 
tion), and interconnections between pores. The number 
of metastable states would however be too large for a 
systematic study. The model was therefore designed to 
exhibit only a few metastable states. 

The main potential sources for generating metastable 
states are nanometer-scale heterogeneities, due to pore- 
size distribution, and variations in surface chemistry. In 
both cases, the prominent effect on the adsorbed fluid 
comes from the modulation of the effective fluid/wall 
interaction. [2^, [2l| We then chose to investigate a cylin- 
drical pore containing few domains of variable fluid/wall 
interaction. The diameter of the nanopore is 8 cr// 
(cT// is the fluid- fluid Lennard- Jones diameter). In this 
way, the smooth-wall approximation can be applied: the 



external potential seen by a fluid particle in the pore is 
calculated by integrating the fluid-wall (6 — 12) Lennard- 
Jones potential (eg/ = 1.277 e//, Ugf = 1.094 ct//) over 
a uniform distribution of substrate interacting sites of 
density 0.8265 ct^j. The calculated reduced external po- 
tential '^*cyi{'r) — '^cyi{r)/eff in this perfectly cylindrical 
pore is given in Fig. [T] The heterogeneity is introduced 
by modulating this external potential along the axial di- 
rection z: 



1 + a(z) sin 



in I 



V L 



Kyiir) (1) 



where a = 0.3 for 2; < and a = 0.2 for z > 0, and 
L = 24(T// is the simulation box length. This external 
potential exhibits four domains of spatial extension of 
a few molecular diameters along the axial direction (see 
Fig. [2]). Periodic boundary conditions are applied along 
the axial direction z. The interactions are truncated at 
half the simulation box size (minimal image convention) . 




FIG. 1: External potential 'I^'j,; (r) for the perfectly cylindrical 
pore as a function of radial distance r. 




FIG. 2: External potential 'iipore{''',z) for the heterogenous 
pore as a function of coordinates r and z. 

The fluid adsorption properties are calculated by 
Monte-Carlo simulations. Thermalization of the ad- 
sorbed fluid is performed by particle displacement trials. 
Chemical equilibration between the adsorbed fluid and 
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the reservoir is done by particle exchange trials. The 
acceptance probabilities are given by the Metropolis al- 
gorithm. The fluid in the reservoir is assumed to be 
ideal, which is a good approximation since the simula- 
tions are performed well below the critical point. As 
a consequence, the gas in the reservoir does not need 
to be treated explicitly, which speeds up the calcula- 
tions. In the limit of infinite reservoir size, one recov- 
ers the usual Grand Canonical Monte Carlo algorithm 
(GCMC). The adsorption/desorption curves are deter- 
mined as follows. The initial configuration is the empty 
pore. Few molecules are then added in one step to the 
system (pore-|-reservoir). In order to mimic the experi- 
mental situation, the molecules are actually added to the 
reservoir cell, which results in an increase in fj.i{, and, 
accordingly, to a lack of equilibrium between the pore 
and its reservoir. The whole system is then relaxed for a 
while. Molecules exchange between the cells, which pro- 
duces a small flow from the reservoir to the pore in the 
first steps. This flow comes to zero when stationarity is 
reached = jip). We emphasize that this matter flow 
does not describe correctly the mass transport, since the 
Monte Carlo algorithm allows transfer of matter every- 
where in the porous substrate. 10^ Monte Carlo trials 
per particle are then performed to acquire statistics for 
computing averaged quantities. This gives the first point 
of the isotherm. The subsequent points are obtained ac- 
cording to the same procedure, i.e. a small increase in 
the number of molecules in the reservoir, followed by a 
relaxation run, and, finally, by a long run for acquisi- 
tion. After complete adsorption has been achieved, the 
desorption isotherm can be calculated by decreasing the 
amount of particles in the reservoir step by step, as in 
a real experiment. This algorithm allows us to obtain 
the scanning curves embedded within the main hystere- 
sis loop. Note that for the case of an infinite reservoir, 
it is the chemical potential which is increased stepwise 
instead of the amount of particles in the reservoir, which 
is then infinite. 



C. Coarse-grained model for adsorption in a 
disordered porous material 

As discussed in previous papers, 0, S l23l | our ap- 
proach to fluid adsorption in disordered mesoporous ma- 
terials is based on a coarse-grained lattice-gas descrip- 
tion which incorporates the essential physical ingredi- 
ents of the solid-fluid system. We consider a three- 
dimensional BCC-latticc of linear size L in which each 
of the Nsp = SL'^ sites may be occupied by a fluid or by 
a solid particle (L is measured in units of the lattice spac- 
ing a and we set a = 1; for a lattice, the volume Vp is then 
simply equal to the number of sites Ns p times the volume 
of the elementary cell = 1, so that we shall speak of 
volume and number of sites indistinctly). Multiple occu- 
pancy of a site is forbidden and only nearest-neighbor in- 
teractions are taken into account. The fluid particles can 



equilibrate, as explained below, with a finite-size reser- 
voir at fixed temperature whereas the solid particles are 
"quenched" and distributed according to a specific choice 
of the porous material structure. For the sake of sim- 
plicity, we study a random matrix with a porosity (i.e. 
fraction of sites without solid particles) = 75%. The 
relevant correlation length of the solid is around one lat- 
tice spacing so that even the smallest systems studied 
(with a linear size L = 25) can correctly describe the col- 
lective effects occurring inside the matrix on long length 
scale (such as a sharp condensation event in the whole 
pore space). 

The starting point of our theoretical analysis is the 
following expression of the free-energy functional in the 
mean-field approximation: 

Fp[{pt}] = ksT^lp., ln/3j + (r/, - /?,,) ln(?7i - pi)] 

i 

- ^ff 5Z P'P^ ~ ^P''^^ ^ + Pi^^ ^ '^')] 

<ij> <ij> 

where pi is the thermally averaged fluid density at site i 
and rji = 0, 1 is a quenched variable describing the occu- 
pation of the lattice by solid particles (rji = if the site i 
is occupied by the solid and rji = I otherwise); Wff and 
Wsf denote the fluid-fluid and fluid-solid attractive in- 
teractions, respectively, and the double summations run 
over all distinct pairs of nearest- neighbor sites. 

We first start with the grand canonical situation {a = 
0) where fluid particles can equilibrate with an infinite 
reservoir that fixes the chemical potential p. Minimiz- 
ing the grand-potential functional rip[{pi}] = Fp[{pi}] ~ 
pJ2i Pi with respect to {pi} at fixed T and p for a given 
realization of the solid yields a set of coupled equations. 



1 + exp[-/3(/x + Wff J2j/, P] + Wsf Ej7.i(l - Vj))] ' 

(3) 

where the sums run over the c = 8 nearest neighbors of 
site i. By using a simple iterative method to solve these 
equations, one finds solutions that are only minima of 
the grand potential surface, i.e. metastable states. For 
a given realization of the solid, the adsorption isotherm 
is obtained by increasing the chemical potential in small 
steps 5 p. At each subsequent p, the converged solution 
a,t p — 6p is used to start the iterations. 

What happens in the mixed situation {a ^ 0)? 
The fluid particles can equilibrate between the solid 
sample and a finite-reservoir so that their total num- 
ber TV is fixed. The isotherm is then obtained by 
increasing in small steps 8N . In this case, the 
system tries to minimize its total Helmholtz free- 
energy FT[{pJ,Pfl,T] = Fp[{p,}, T] + FflK, T], where 
Fp is the Helmholtz free-energy of the gas reservoir 
{Fr\pr,T\ = VR{kBT\pR\iipR + (1 - pfl)ln(l - pr)\ - 
Wffcp^R/i} with Vr the volume of the reservoir, i.e., 
again up to a^^ = 1, the number of sites of the reservoir, 
while satisfying the global constraint A^ = Np + Nr ~ 
X]j Pi + Vrpr. This can be solved in a natural way by 
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the method of Lagrange multipUers. We consider the III. RESULTS ON RESERVOIR-SIZE 

function DEPENDENCE FOR THE ATOMISTIC MODEL 



nT[{pi\, PR, X, T] = Ft[{p,}, PR, T]+\{N~Y.p,-Vrpr}, 

i 

(4) 

where A is a Lagrange multipUer that has the meaning of 
the chemical potential coupled to the densities. Minimiz- 
ing Ft with the constraint on densities amounts to solve 
simultaneously the coupled equations = 0, = 

and = or equivalently 

kBT\n[ ^' ] - A - Wff^pj - Wsf ^(1 - ijj) = 0, 
iji Pi ... ... 

fcBTln[ ] - A - WffPR = 0, 

- PR 

N -Y,P,~Vrpr = Q, l<i<Nsp. (5) 

i 

One has then to define an iterative scheme that specifies 
how the system goes from one converged solution to an- 
other as the total number of particles is slowly changed. 
The details are given in the appendix. 

We are searching for configurations for the local den- 
sities Pi that are local extrema of the grand-potential 
functional ilp for the special value of the chemical po- 
tential fi ~ X. However, it is not fully guaranteed that 
the above algorithm necessarily converges to a local min- 
imum, nor even to an extremum, since the constraint 
could in principle stabilize an unstable state. Therefore, 
we regularly ascertained that configurations obtained in 
the mixed ensemble are stable in the grand-canonical en- 
semble, i.e. that they indeed are metastable states, by 
starting grand-canonical calculations with the converged 
solution at chemical potential A. In addition, it must be 
emphasized that the global constraint is not satisfied un- 
til convergence is reached, in the spirit of the Lagrange 
method. It is therefore doubtful that one can attribute 
any physical meaning to the intermediate stages of the 
iterative process. Recall also that one does not take into 
account mass transport: local densities can change ev- 
erywhere and instantaneously inside the porous sample. 

Moreover, we focus on the adsorption (fiUing) process. 
As shown in previous papers, [sl. |24|| fiuid desorption may 
crucially depend on the presence of an external surface 
for the porous solid: the system then includes a real inter- 
face between the adsorbed fluid and the external vapor, 
and. during desorption, the vapor domain may penetrate 
and drain the solid from the outer surface (the so-called 
percolation and depinning transitions discussed in Ref. 
[24| . These mechanisms are no longer bulk phenomenon 
and their, a priori more subtle, analysis will be carried 
out elsewhere. We have therefore used periodic boundary 
conditions for the sample and the reservoir (separately). 
The ratio y = Wgf/wff has been fixed to 0.8 so that 
the adsorption isotherms exhibit the interesting range of 
phenomena as the temperature changes (see Ref.Q). 



A. Grand-canonical isotherms 

We first focus on the adsorption/desorption isotherms 
obtained in the limit of infinite reservoir size (a = 0), i.e. 
the GCMC data. The results are given in Fig. [3] for var- 
ious temperatures. As can be seen, the highest temper- 
ature isotherm (reduced temperature T* = ksT/eff = 
1.0) is perfectly reversible, i.e. the adsorption and des- 
orption curves are superimposed. On the other hand, for 
lower temperatures, the adsorption and desorption curves 
differ and exhibit hysteresis with vertical steps. We em- 
phasize that the vertical lines are guides to the eye and 
do not correspond to GCMC data. The isotherm jumps 
in one step, and cannot be stabilized in between. 

Quite noticeably, the simulation points may be 
grouped into branches. For instance, for a reduced tem- 
perature of 0.9, three branches are present: the gas-like 
branch (with a fiuid layer adsorbed at the wall, the rest of 
the pore being filled with gas), which exists down to very 
low /i, the liquid-like branch (pore filled by dense fluid) 
at high /i, and a branch of intermediate density which 
exists only for a limited range of p. We have found that 
adsorption and desorption are reversible along any given 
branch. These branches are associated to deep local min- 
ima in the free energy landscape describing the system 
(the grand potential). The energetic barriers separating 
these minima arc generally large compared to the ther- 
mal fluctuations sampled by the Monte-Carlo algorithm: 
the system then remains trapped in these local minima, 
which explains the reversibility of adsorption/desorption 
along the branches. However, for some particular values 
of p the barriers become sufficiently small for allowing 
thermal fluctuations to make the system jump into an 
adjacent local minimum. These values dcflne the limit 
of stability of the branches (in the grand canonical en- 
semble). The vertical lines in Fig. [3] indicate the new 
local minimum (branch) reached by the system. Note 
that the ^-range of the various branches overlap, which 
means that for some p, the system may be stabilized in 
GCMC at various degrees of pore flUing p. This is the 
main origin of hysteresis. As can be seen, the lower the 
temperature, the larger the /x-range of existence of the 
branches, and the wider the hysteresis. The number of 
branches increases for decreasing temperature, at least 
down to T* = 0.8 (Fig. [S]). We shall show later that this 
is actually still true down to T* = 0.6 if one takes into 
account metastable states that are hidden in the GCMC. 
For T* = 0.9, the three branches are visited during both 
adsorption and desorption. They belong to the main 
adsorption and desorption curves. For T* = 0.8, flve 
metastable states belong to the main adsorption curve, 
and three to the main desorption curve. The branch of 
intermediate density reached during desorption (second 
branch, p* = pcr^f ~ 0.4) looks like a continuation of 
the third branch visited during adsorption (p* ~ 0.4). 
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This point was checked by showing the reversibihty of 
the complete branch by GCMC (by increasing and de- 
creasing fj,). The total number of metastable states is 
then five, and their range of existence was determined 
in a systematic way by /i-variations. The results are 
shown in Fig. 21 Here again, lines are guides to the 
eye, and the vertical lines show the new metastable state 
reached by the system beyond the stability limit. The 
case T* = 0.6 deserves special attention. The procedure 
that consists in following simple ascending and descend- 
ing /x-paths does not allow to reach the two branches 
with reduced densities around 0.3 and 0.45. These states 
were obtained from the corresponding states at higher 
temperature T* = 0.8 by following a particular (/-i,T) 
path along which their stability is preserved. This path 
essentially consists of slowly decreasing the temperature 
and increasing the chemical potential in order to keep 
the average number of particles roughly constant. When 
T* = 0.6 has been reached, the /i-range stability can be 
determined for both states (see Fig. [J]). It is possible 
that more metastable states actually exist, but the pro- 
cedures previously explained did not allow finding more 
than five branches in our simple system. 




FIG. 3: Adsorption/desorption isotherms obtained by GCMC 
simulation for various temperatures (symbols). Lines are 
guides to the eye. 



B. Metastable states 

Figure [5] displays an enlargement of the lowest tem- 
perature (T* = 0.6) GCMC results, without the vertical 
lines. As previously mentioned, the simulation points 
group into reversible branches (solid lines are guides to 
the eye) which correspond to local minima in the free- 
energy landscape. We now give a molecular-level descrip- 
tion of these local minima. The lowest density branch 
corresponds to gas-like fluid filling the pore with ad- 
sorbed molecules at the wall. The highest density branch 
corresponds to liquidlike fluid saturating the pore. The 
three other states are stabilized by the chemical corruga- 




FIG. 4: Stables and metastables states obtained by GCMC 
simulation for various temperatures (symbols). Lines are 
guides to the eye indicating the path followed by the system 
beyond the stability limits of each metastable state. 

tion. Visual inspection of molecular configurations (see 
lower panel of Fig. [5]) shows that the three branches 
correspond respectively to: one liquid-like meniscus in 
the most attractive domain of the pore; two liquid-like 
menisci in the two attractive regions of the pore; one sin- 
gle large meniscus bridging the two attractive domains. 



C. Mixed-ensemble results 

We now focus on the adsorption/desorption results ob- 
tained for a finite-size reservoir. Two ratios Vj^/Vp = 500 
and 2000 have been considered (a ~ 2. x lO^'^ and 
5. X 10~^). The procedure previously used to obtain all 
metastable states will be applied here in its simplest form, 
i.e. isothermal /i-paths, because it corresponds to what is 
actually feasible for realistic mesoporous substrates and 
in real experiments (one essentially measures the main 
adsorption/desorption hysteresis and scanning curves). 
The results for the main hysteresis loops are given in 
Fig. [SI where the grand canonical case {a = 0) has also 
been plotted for comparison. In all cases the symbols are 
the simulation data, and the color (thick) lines are guides 
to the eye that connect the simulation points in the same 
order as they were obtained by slowly increasing (ad- 
sorption) or decreasing (desorption) the total number of 
particles in the system. The thin (black) lines are the five 
branches previously obtained by extensive GCMC study 
(Fig. [5]). In the limit of infinite reservoir size {a = 0), 
starting from the empty system and gradually increas- 
ing fi results in the continuous filling of the system up 
to fi* = fi/cf f = —9.43 where the system jumps into the 
second local minimum. Increasing further the chemical 
potential makes the system jump directly to the fifth lo- 
cal minimum. Upon desorption, the system remains in 
the saturated state for a while (hysteresis) and finally 
jumps (for fi* = —9.94) directly into the first local min- 



7 



? 0.4 - 





FIG. 5: Enlargement of the low-temperature metastable 
states found by GCMC. Solid lines indicate reversible 
branches. The lower panel shows at ^* = —9.52 (dotted line), 
for each of the five metastable states sketched in the upper 
panel, one typical molecular configuration. 



imum (gas-like branch) without visiting any intermedi- 
ate state. The vertical lines represent the constraint im- 
posed by the reservoir (/x*= cte) but not necessarily the 
path actually followed by the system during the transi- 
tion (transient phenomenon are not properly described 
by the Monte-Carlo algorithm). 



In the cases of finite reservoir, it is the total number 
of particles in the system plus its reservoir which is in- 
creased by small steps as in a real experiment. As can 
be seen from Fig. [HI the amount of adsorbed fluid ini- 
tially increases along the gas-like branch in a continous 
manner until it reaches its stability limit. (Note that 
this limit occurs at slightly larger ^ for smaller reservoir: 
-9.43 for a = 0, = -9.41 for a ~ 5. x lO'^ and 
^* = -9.36 for a ~ 2. X 10"^.) At this stability limit, a 
small addition of extra molecules in the reservoir destabi- 
lizes the adsorbed fluid which then jumps into the second 
state. During the transition, particles leave the reservoir 
and adsorb into the pore in the most attractive region 
to form a meniscus (second metastable state), which re- 
sults in a decrease of the chemical potential. Here again, 
our algorithm is not meant to describe this transition. 
The path actually followed by the system is not known. 
The (curved) line shown in the flgure corresponds to the 
constraint of conservation of the total number of parti- 
cles in the system plus its reservoir. After equilibration, 
the system reaches the second branch. Further increase 
in the total number of particles makes the system follow 
this branch until it reaches a new stability limit. (As pre- 
viously, this limit slightly increases for a finite reservoir: 
/I* = -9.27 for a = 0, //* = -9.23 for a ~ 5. x 10~^ 
and /i* = -9.24 for a ~ 2. x 10"^.) The system then 
jumps onto the fifth branch for a = 0, onto the fourth 
branch for a ~ 5. x 10^"*, and onto the third branch for 
a ~ 2. X 10-3. 

The above results illustrate one of the main effects of 
the reservoir: the system does not necessarily visit the 
same metastable states as one changes the size of its 
reservoir. From Fig. [6l this effect may be tentatively 
interpreted as follows: the reservoir imposes a constraint 
on the relation between the chemical potential and the 
amount adsorbed, materialized by the vertical or inclined 
curves. After a branch stability limit has been reached, 
the system approximately follows this line of constraint 
during the transition, and finally meets a new branch 
at a location that depends on the ratio between the 
reservoir and the system size. To illustrate this point, 
starting from the highest /i GCMC point of the second 
branch (limit of stability), we have determined the new 
metastable state reached by the system after a small in- 
crease in the total number of molecules for various reser- 
voir sizes. The results are given in Fig. [71 This picture is 
however somewhat oversimplified and does not apply in 
all situations since it does not take into account the com- 
plexity of the underlying free-energy (grand-potential) 
landscape that determines which path is actually followed 
by the system (being in configurational space, the free 
energy landscape is an object of very high dimension). 
At a molecular level, this means that complex nucleation 
and fluid adsorption processes determine the path actu- 
ally followed by the system. One of the consequences 
is that the system does not necessarily stop on the first 
branch crossed by the constraint line. This is illustrated 
in our simple system for a ~ 2. x lO"'^ during desorp- 



8 



tion from the fourth branch: as can be seen, the con- 
straint Hne crosses the third metastable branch, but the 
system avoids this state and actually reaches the second 
branch. In the grand-potential landscape picture, this 
means that, starting from the fourth local minimum, and 
following the steepest slope, the system has not crossed 
the basin of attraction associated to the third minimum 
and has followed its way until being finally trapped by 
the basin associated to the second local minimum. In 
the molecular-level description, the fourth branch corre- 
sponds to a large meniscus bridging the two attractive 
regions. Upon desorption, the system would have to nu- 
cleate a bubble in between the attractive regions in order 
to reach the third branch corresponding to two menisci 
in the attractive regions (see Fig. [8l upper path). This 
nuclcation barrier is probably too high, and it is more 
favorable for the system to dcsorb by recession of the 
meniscus in the less attractive regions until it is left with 
one single meniscus in the most attractive region corre- 
sponding to the second state (see Fig. [51 lower path). 

As for the grand-canonical case, it is possible to per- 
form a systematic search for all possible metastable states 
in the system in contact with a finite reservoir. The five 
branches already described are found in all cases. Fig- 
ure [9] shows the superimposition of all simulation points 
(stable over long runs) obtained for the various ratios be- 
tween the reservoir and system sizes already presented. 
As can be seen, the points group onto branches corre- 
sponding to the five metastable states found in the sys- 
tem. It is however reminded that the limits of stability 
of the branches slightly vary with the reservoir size. This 
may be related to the fact that transitions proceed via 
energetic barriers and that the amplitude of the fluctua- 
tions allowed by the reservoir depends on its size. To 
summarize the results, reducing the size of the reser- 
voir allows one to explore metastable states upon ad- 
sorption/dcsorption that are not accessible in the grand- 
canonical situation. It is also found that the number of 
metastable states visited by the fluid increases with de- 
creasing reservoir size. 



IV. RESULTS FOR THE COARSE-GRAINED 
MODEL 

A. Grand-canonical isotherms 

Typical results for adsorption isotherms obtained at 
various temperatures in the grand-canonical ensemble are 
shown in Fig. 1101 When the temperature decreases, 
the isotherm shape changes from smooth to steep. As 
discussed in previous papers, 2J] this corresponds to a 
true out-of-equilibrium phase transition, the so-called 
"avalanche transition" ,[7[ with the sudden appearance 
of a macroscopic, connected liquid domain in the whole 
porous sample. A previous extensive scaling study [24 1 
showed that it exists a critical value of the temperature 
for which the isotherm changes from continuous to dis- 




FIG. 6: Adsorption/desorption isotherms for various values 
of the relative size of the gas reservoir: q = 0, a ~ 5. x 10~* 
and 2. x 10"^, where a = Vp/{Vp + Vr). 



continuous in the thermodynamic limit. At higher tem- 
perature, isotherms look gradual but at reduced temper- 
atures T* = 1.40, they consist of little steps of varying 
sizes (see insert in Fig. [TT|) . 

The above results are related to the characteristics of 
the grand-potential landscape. At low temperatures, this 
multidimensional landscape (recall that it is a function 
of the local fluid densities, with here Nsp ~ 10^ — 10^) 
is characterized by a large number of local minima, the 
metastable states in which the system may be trapped. 
Since thermally activated processes are neglected in the 
LMFT (an approximation that, as already discussed, 
finds its justification in the experimental reproducibil- 
ity of the adsorption isotherms on the time scale of most 
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FIG. 7: Circles: GCMC results for T* = 0.60. Lines are 
guides to the eye. The arrows indicate the jump made by the 
representative point of the adsorbed fluid for various reser- 
voir sizes, after destabilization of the initial point A by a 
small increase in the total number of particles in the system. 
The curved arrow indicates increasing values of the relative 
reservoir size: a = 0. (GCMC), 2.2 x 10"*, 2.5 x 10"", 5.0 x 
10"", 1.0 X 10"^ 2.0 X 10"^ 2.4 X 10"^. 
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FIG. 8: Possible mechanisms for the evolution of the fluid 
conflguration inside the pore starting from a given metastable 
branch (see text). 



experiments), the evolution of the system is only due to 
a variation of the external driving (here the chemical po- 
tential). As /X varies, the system either follows the min- 
imum in which it is trapped as this minimum deforms 
gradually (the flat portions in the insert of Fig. [TT|) . or 
it falls instantaneously into another minimum when the 
former reaches its stability limit. This later move is a dis- 



FIG. 9: Metastable states of the adsorbed fluid obtained for 
all the different sizes of the reservoir studied here (shown with 
different colors; colors on line). 
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FIG. 10: Hysteresis loops around the "capillary condensa- 
tion" in the grand canonical ensemble for a L = 100 sample 
as a function of the temperature: T* = 1.40 (full blue line), 
T* = 1.10 (dotted red line), T* = 0.80 (dashed black line) 
(colors on line). For the last two temperatures, only part of 
the desorption branch is shown. 



continuous and irreversible process, an avalanche, which 
is at the origin of the history-dependent behavior of the 
system, e.g. the hysteresis. The avalanche corresponds 
to some collective condensation event inside the porous 
sample that manifests itself by a jump in the adsorption 
isotherm. The size of the avalanche may be macroscopic 
as for T* = 0.8, and the adsorption isotherm is discon- 
tinuous in the thermodynamic limit; or it may be mi- 
croscopic as for T* = 1.40, and the adsorption isotherm 
remains smooth in the thermodynamic limit. (More cal- 
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culations should be necessary to conclude on what hap- 
pens at T* = 1.10.) The hysteresis loop encloses all 
the metastable states of the system: this is illustrated in 
Fig. [TT] by the scanning curves obtained by performing 
incomplete filling of the matrix and then decreasing the 
chemical potential to drain the adsorbed fluid. Of course, 
the number of metastable states is very large, probably 
exponential in the system size, and only few of them are 
revealed with this simple procedure. 




FIG. 11: Grand canonical descending scanning curves in a 
L — 50 sample at T* = 1.40. A close-up of the adsorption 
isotherm is shown in the inset. 



B. Canonical and mixed-ensemble isotherms 
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FIG. 12: Canonical (symbols) and grand-canonical (lines) 
adsorption isotherms in a L = 100 sample at T* = 1.10 (right) 
and T* = 0.80 (left). 

What happens when the porous sample is coupled to a 
finite reservoir? First, we begin with a vanishingly small 



reservoir, which is the canonical situation where one con- 
trols the number of adsorbed particles. The main results 
arc summarized in Fig. [T2] where we compare the ad- 
sorption isotherms obtained with the "/i-driven" (grand- 
canonical) and "p-driven" (canonical) procedures. The 
behavior is quite different in the low and the high temper- 
ature regimes, respectively characterized by the absence 
and the presence of a macroscopic /x-driven avalanche. 
At the highest temperature studied (T* = 1.80, not 
shown in the figure), there exists only one stable state, 
the equilibrium one, and controlling either the adsorbed 
density or the chemical potential yields the same result. 
This is also true for all temperatures at very low ad- 
sorbed densities (not shown in the figure). At lower 
temperatures, avalanches appear in the grand-canonical 
isotherm. Collective localized events, which we define 
as "avalanches" irrespective of the control variable, also 
show up in the canonical isotherm in the form of small 
jumps in the chemical potential toward a lower value: 
with the disappearance of the initial minimum, the sys- 
tem has to find another metastable state with the re- 
quired adsorbed density; as this state cannot be found 
at a higher value of the chemical potential since the 
grand-canonical adsorption/desorption isotherms have 
been shown to represent the extremal curves that en- 
compass all the metastable states of the system,]^ the 
chemical potential must decrease. (Note that in a grand- 
canonical setting "avalanches" appear as jumps in the ad- 
sorbed density p at constant chemical potential fi whereas 
in a canonical one, they appear as jumps in fi at constant 
p.) When a metastable state is found and fluid is further 
added, the system smoothly follows this state (see the 
quasi-linear portions in Fig. [T3| until it reaches the cor- 
responding stability limit. Then, there is a new jump in 
the chemical potential to a smaller value, and the evo- 
lution proceeds in this way until the porous sample is 
completely filled with liquid. The contrast between p- 
driven and p-driven isotherms that is illustrated in Fig. 
1121 is very reminiscent of what was found in ferromag- 
netic systems when comparing rnagnetization-driven and 
magnetic-field-driven protocols. |13| As in the latter case, 
the canonical (/9-driven) isotherms are closely related to 
the distribution of metastable states inside the (grand- 
canonical) hysteresis loops. 

We now discuss in more detail the canonical (p-driven) 
isotherms. At T* = 1.40 with a sample of linear size 
L = 100, the jumps in p(p) are very small and one 
needs to zoom in on the isotherm to see small differ- 
ences between the canonical and grand-canonical proto- 
cols, as illustrated in Fig. [T31 This smallness is related 
to the presence of many metastable states in the close 
vicinity of the grand-canonical isotherm. Since both the 
grand-canonical and the canonical avalanches remain of 
microscopic size, we expect that the small jumps in ei- 
ther p or p become infinitesimally small in the thermo- 
dynamic limit. This is indeed what is observed in Fig. 
1131 where isotherms obtained for different system sizes 
are compared. Therefore, we predict that the canonical 
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FIG. 13: A close-up of the canonical (symbols) and grand- 
canonical (lines) adsorption isotherms at T* = 1.40 for differ- 
ent linear sizes of the system (colors on line). 



and grand-canonical curves should become identical in 
the thermodynamic limit. 

A similar behavior is observed at T* — 1.10, as il- 
lustrated in Fig. 1141 However, the jumps in p are larger 
than at T* = 1.40 and are clearly associated with a much 
more fluctuating chemical potential each jump in 

p generates a re-entrant behavior in p{p). It is difficult 
to conclude on how the isotherm will evolve in the ther- 
modynamic limit: as the size of the system increases (see 
Fig. I14p , there are more and more intermediate points in 
the steepest part of the grand-canonical isotherm and. at 
the same time, the canonical isotherm is less and less fluc- 
tuating. The temperature T* ~ 1.10 is probably barely 
above the critical temperature of the (grand-canonical) 
avalanche transition and we expect that both the canoni- 
cal and the grand-canonical isotherms would become con- 
tinuous and would coincide in the thermodynamic limit. 

On the other hand, T* = 0.80 is undoubtedly below 
the critical temperature of the avalanche transition: all 
samples studied, whatever their size, display a large jump 
in their grand-canonical isotherms (see Fig. [T5|) . (In 
a previous work on a random matrix. [23| we concluded 
through an extensive flnite-size study that a macroscopic 
jump exists in the thermodynamic limit at T* = 0.8 
for y = 0.9; since the value y = 0.8 used in this pa- 
per corresponds to a weaker disorder than y = 0.9, this 
is necessarily true here as well.) As shown in Fig. I15[ 
the canonical curves show a pronounced reentrant be- 
havior, with a large difference with the grand-canonical 
ones: a whole region void of metastable states then ap- 
pears. There are always finite-size effects, which affect 
both the canonical and the grand-canonical isotherms. 
In the latter case, the rare events that trigger the macro- 
scopic avalanches are very sensitive on details about the 



structure of the matrix, and the corresponding chemical 
potentials vary strongly with the system size at small 
L. We have nonetheless checked by performing grand- 
canonical isotherms on very large systems of linear size 
L = 200 that the position of the jump in the thermo- 
dynamic limit is very close to the jump found in the 
isotherm for L = 100 that is displayed in Fig. [T21 On 
the other hand, even if the chemical potentials of the 
canonical isotherms fluctuate less and less as the size of 
the system grows (especially at the end of the adsorption 
process, when the porous sample is nearly filled with liq- 
uid), the overall location of the isotherm in the p-p plane 
does not change significantly with L. 

To discuss the fluctuations of the chemical potential 
with the density in the canonical protocol, we have per- 
formed flnite-size studies of the flrst moments of their 
distribution, after averaging over disorder (i.e. porous 
sample) realizations. More precisely, we compute the 

mean p{p) (hereafter an overbar denotes the average over 

2 

disorder) and the variance A^(p) = p{p)'^ — p{p) for 
three sizes of the system. As explained in Rcf. [l^ . 
the standard argument concerning self-averaging quan- 
tities, whose value in a macroscopic sample is equal to 
the average over all disorder realizations, cannot be ap- 
plied: as a consequence, one cannot be sure that p{p) 
is a self-averaging quantity. To delve more into its be- 
havior, it is interesting to investigate the dependence on 
system size of the standard deviation cr^(p) = A^(p). 
The result is shown in Fig. [TH where (t^ (/9) is plotted as 
a function of the adsorbed density (computed in bins of 
width 0.01). Except at low and high density, the stan- 
dard deviation remains nearly constant with p (typically, 
between p = 0.2 and p = 0.6). It decreases with L but 
quite slowly: we flnd that in the region between p = 0.2 
and p = 0.6; A^(L) = a^{LY ^ L~'' with a flnite-size 
scaling exponent 7 around 1.4, i.e. signiflcantly smaller 
than 3. This indicates that u{p) is self-averaging, but 
only weakly so. As in Ref. [I4I, we have checked that the 
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FIG. 14: Canonical (symbols) and grand-canonical (lines) ad- 
sorption isotherms at T* = 1.10 for samples of different sizes: 
L = 25 (left), L = 50 (middle), L = 100 (right). 
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FIG. 15: Canonical (symbols) and grand-canonical (lines) ad- 
sorption isotherms at T* = 0.80 for samples of different sizes: 
L = 25 (left), L = 50 (middle), L = 100 (right). 




FIG. 16: Left: Standard deviations of the chemical poten- 
tials as a function of the adsorbed density p at T* — 0.80 for 
different system sizes: L — 25 (black), -L = 50 (red), L = 100 
(blue). Lines show averaged values. Right: histograms of 
fi{p) for different values of p (p = 0.2 in red/left, p — 0.6 in 
black/right ) and different system sizes {L = 25: top, L — 50: 
middle, L — 100: bottom). 



corresponding histograms are roughly Gaussian (see Fig. 
[16]). Although the situation with respect to self-averaging 
remains somewhat unclear as far as the whole isotherm 
is concerned. [l3| it is nonetheless instructive to study the 
mean fj,(p) as a function of system size. This is displayed 
in Fig. I17[ where one can see that the /x(p)-isothcrms 
vary significantly with the size of the system, becoming 
steeper as L increases. We cannot conclude whether or 
not the isotherm becomes vertical in the thermodynamic 
limit, but the data strongly suggest that it will be dis- 
tinct in this limit from the grand-canonical isotherm: the 
/i(/9)-isotherms for all studied system sizes intersect near 
the point fi* — —4.47 and p = 0.3, at a significantly 
lower chemical potential than those found for the grand- 
canonical jumps {p* ~ -4.43 - 4.44) for L = 200 (see 

Fig. [IID. 

We now comment on the two types of behavior found 
along the p-driven adsorption isotherm. For increasing 
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FIG. 17: Average canonical adsorption isotherms p(p) at 
T* — 0.80 for different system sizes; from left to right in 
the upper part of the plot (colors on line): L = 25 (black), 
L = 50 (red), L — 100 (blue). The corresponding numbers 
of realizations are 256, 64 and 8, respectively. Typical grand- 
canonical adsorption isotherms for L = 200 samples are also 
shown (dashed curves). 



adsorbed density p, the chemical potential may cither in- 
crease continuously or decrease by discontinuous jumps. 
In the former case, there are only slight swellings of the 
liquid domain, with the liquid-gas interfaces retaining the 
same shapes locally and all the p.^'s increasing. In the 
latter case, we find that each jump corresponds to a sin- 
gle condensation event, a p-driven "avalanche", in which 
all sites of a compact region become liquid. This is il- 
lustrated in Fig. [iHl Considering the two consecutive 
configurations with average densities p and p -I- Ap (be- 
fore and after the jump), sites are considered as turning 
liquid (respectively, gas) when the variation of the local 
fluid density is larger than 0.3 (respectively, lower than 
—0.3) and are shown in red (respectively, grey) in the 
figure. The size of the condensed region varies for each 
jump along the isotherm, but the number of particles in- 
volved in the local condensation can go much beyond the 
controlled increment VpAp: in Fig. [181 the condensation 
event (in red) corresponds to a local increase in the num- 
ber of particles that is 110 times larger than the initial 
increment. Therefore, others regions of the porous sam- 
ple must be drained: as shown in the figure, the most 
important evaporation events (in grey) remain smaller 
and less compact than the condensation event and ap- 
pear at random throughout the material. (Whereas the 
definition of "turning liquid" is well characterized since 
the histogram of positive variations of the local density 
is always peaked at a value close to 1, this is much less so 
for its counterpart "turning gas" : the histogram of ncga- 
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FIG. 18: Condensation and evaporation events between p = 
0.3962 and p — 0.3963 in the p-driven protocol. 

tive variations is continuously decreasing between and 
— 1, so that the decrease of the average density that bal- 
ances the condensation event mostly conies from small 
but widespread local variations.) This suggests the fol- 
lowing interpretation: adding a small amount of fluid 
in the sample may trigger a condensation "avalanche" 
whose size is larger than the added fluid amount; a de- 
crease of the chemical potential allows a slight draining of 
the rest of the sample and provides the required amount 
of fluid to "feed" the avalanche. 

Finally, Fig. \W\ presents our results in the mixed- 
ensemble for different sizes of the reservoir. As the size 
increases, the /x(p)-isotherms exhibit flat segments with 
decreasing length and jumps in the fi-p plane that arc 
more and more pronounced. This behavior is similar, al- 
beit in a more complex systems and in the presence of a 
much larger number of metastable states, to that found in 
the atomistic model and illustrated in Fig. [S] When the 
volume of the reservoir becomes 2000 times greater than 
the volume of the porous material, the isotherm displays 
a unique jump between the grand-canonical low-density 
branch and the grand-canonical high-density branch, but 
it still comes with a decrease of the chemical potential 
(lower right panel of Fig. [T9|) . For larger reservoir sizes, 
the jump gets closer and closer to the grand-canonical 
isotherm. Convergence however is rather slow. On the 
other hand, for the smallest reservoir sizes (a reservoir 
volume 500 times the volume of the porous material or 
less: see the two upper panels of Fig. [T9|l . the isotherms 
tend to superimpose on the canonical isotherm. It can 
be however noticed that the explored metastable states 
are not necessarily identical. As explained in section III, 
one indeed has to be very careful with the representation 
of states in the fi-p plane which is only a projection from 



the highly multi-dimensional configurational space: it is 
possible to cross branches as the mixed-ensemble con- 
straint allows the system to explore other areas of the 
underlying free-energy landscape. 




FIG. 19: Adsorption isotherms at T* = 0.80 for different 
relative sizes of the reservoir for a L = 50 sample. In the 
four panels, the canonical isotherm is shown in light grey and 
the grand-canonical isotherm as a dashed curve. Top left: 
a = 10"^, top right: a = 2. x 10"^, bottom left: a = 10"^, 
bottom right: a = 5. x 10^"', where a = Vp/{Vp + Vr). Note 
that both the slope and the length of the inclined segments 
that represent jumps in the p-p plane (the lines are therefore a 
guide to the eye) increase as a decreases: the slope is directly 
related to a (compare with Fig. [Tjl. 



V. CONCLUSION 

In this paper, we have presented a study of the out- 
of-equilibrium, hysteretic response of a fluid adsorbed in 
inhomogeneous porous materials when one couples the 
porous sample with a flnite-size reservoir and controls 
the total number of particles. Varying the relative size 
of the sample and the reservoir allows one to interpolate 
between a canonical situation with a controled adsorbed 
density and a grand-canonical situation with a controled 
chemical potential. We have considered both an atom- 
istic model of a fluid in simple, yet structured pore and a 
coarse-grained model for adsorption in a disordered meso- 
porous material. The adsorption isotherms have been 
computed by molecular Monte Carlo simulation in the 
former case and by a density functional approach in a 
local mean-field approximation in the latter. The Monte 
Carlo simulations give us a clear picture of what occurs 
at the molecular scale in a small and weakly disordered 
system whereas the density functional approach provides 
insights at the mesoscopic scale for a fully disordered sys- 
tem. In both cases, we have found metastable states that 
appear as branches of finite extent in the p^p plane. The 
metastable states correspond to inhomogeneous configu- 
rations of the fiuid and the number of branches increases 
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rapidly with the complexity of the material. It is worth 
stressing two points: first, these states are not unsta- 
ble, i.e. not purely stabilized by the constraint on the 
total number of particles [25j: second, metastability is in- 
duced by the intrinsic inhomogeneity of the solid. Ac- 
cordingly, the picture of metastability in such systems is 
quite different than that encountered in mean-field-like 
descriptions of homogeneous (bulk) systems undergoing 
first-order transitions. 

We have shown that the way the system evolves be- 
tween these metastable states may depend on the pro- 
tocol, controlled here by the relative size of the reser- 
voir. In particular, our results suggest that a discon- 
tinuity in the grand-canonical adsorption isotherm (an 
out-of-equilibrium "avalanche transition") is associated 
with the absence of metastable states in a whole region 
of the /i — jO plane and that the corresponding canonical 
adsorption isotherm (and, more generally, isotherms per- 
formed with a small enough reservoir size) differs from 
the grand-canonical one and displays a reentrance, even 
in the thermodynamic limit. 

What is the relevance of our results for experimen- 
tal set-ups ? For illustration, we consider two examples. 
In order to determine the adsorption isotherm with the 
volumetric method using for instance nitrogen, known 
amounts of fluid are admitted stepwise in the sample cell. 
The amount of adsorbed fluid is the difference between 
the admitted gas and the amount of gas that fills the 
"dead volume" , i. e. the free space in the sample cell in- 
cluding connections: this is equivalent to the procedure 
discussed in this paper, with the dead volume playing the 
role of a finite-size reservoir for the sample. The amount 
of gas in this reservoir is calculated from the fiuid equa- 
tion of state and from measurements of the pressure, the 
temperature and the dead volume. In an experiment with 
a sample size of 20 mm"^, a cell of a few cm^ could be used. 
In such a situation, the ratio of volumes (porous medium 
versus total) is around 10"^ and the ratio of adsorbed 
amount on the total amount around 10^^. On the other 
hand, to measure Helium adsorption in aerogels, E. Wolf 
et al [1^ use a different experimental set-up that was 
first proposed by Herman et alQ, in which the adsorbed 



amount is controlled through the temperature of a helium 
gas reservoir connected to the experimental cell. The to- 
tal amount of Helium is fixed, but varying the temper- 
ature of the reservoir transfers atoms from the reservoir 
to the cell, or conversely. In the recent experiments . [Toj 
with temperatures ranging from 4K to 5K, the adsorbed 
amount at saturation corresponds to around 20% of the 
total amount. 

In both above examples, the ratio between the ad- 
sorbed amount and the total amount of fluid is fully 
in the range of the parameters of our study. However, 
in most experiments in disordered porous materials the 
adsorption isotherms display hysteresis but are smooth 
and continuous: in consequence, as predicted by the 
present study, no effects of the reservoir size are ex- 
pected for macroscopic systems. We have shown here 
that for such effects to be observable, the temperature 
should be low enough for the grand-canonical adsorption 
isotherm to exhibit a true discontinuity and not only a 
very steep variation. It is unlikely that this true discon- 
tinuity could be observed in common disordered porous 
materials such as Vycor or xerogels: the disorder and 
the confinement are too strong and prevent the appear- 
ance of the avalanche transition at temperatures higher 
than the triple point. The case of helium adsorption in 
aerogels is however more promising since at low temper- 
ature (below 4K), jumps have been predicted p?!. [28jand 
may have already been observed in some experiment Q. 
It would therefore be interesting in this case to perform 
controlled experiments for different relative sizes of the 
reservoir to see if the predicted reentrance of the isotherm 
as one moves away from the grand-canonical situation is 
indeed encountered, as observed in hysteretic martensitic 
transformations. 11 
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APPENDIX: DFT ALGORITHMS 



As a starting point, it is convenient to rewrite the con- 
dition of minimization of the grand potential, — g 
as 



exp(-/3A)pi = (r?,-p,)exp[/3(«;//^Pj-|-w^/^(l-?7j))] 

(A.l) 

and sum over i so as to express A as a function of the 
densities 

exp(-/3A) = 

E - Pi) exp[/3(u;// Ej/,; Pj + Wsf Ej^ll - %))]} 



E» P^ 



(A.2) 



Our algorithm is then the following: changing the total 
density px by a small step Ap, i.e. p^^™ = p^'' -I- Ap, 
one first supposes that the supplementary amount is ad- 
sorbed in the porous material: p^^"^ = p"''' -|- Ap. One 
then computes the new Lagrange parameter from 



exp(-/3A"""') = 

Ea('y. - pf) cxp[/3K/ pf + ^sf - %■))]} 



Vpp" 



(A.3) 



(with the old local densities) and the new local densities 
from 



1 + exp[-/3(A-- + wff Pf + - Vi))] ' 

(A.4) 

The new reservoir density p^^ is obtained by inverting 



„new 

Pr 



1 



Pr 



U'ffPR 



(A.5) 



One then checks if the constraint is satisfied, i. e. if pr = 
a ^'y' 1- (1 — oi)p^R^ within a given precision and, if 



not, one iterates using Ap = ^[pi 



Vp 



(1 



a)p^^^'] until convergence is reached. 

In practice, to improve the convergence of the iteration 
procedure, we use a mixing scheme, retaining a part of 
the previous iteration for the subsequent iteration. It 
appears that the configurations visited do not depend on 
the mixing parameter when the convergence criterions 
are strong enough. 

Note that other equivalent algorithms can be devised 
as well and, interestingly, the output of the calculation 
appears to be quite robust to the choice of the algorithm. 
For instance, changing Eq. IA.2I by the equivalent equa- 
tion. 



exp(— /3A) 



Y.i{v^- P^) 



Y.^{P^ CXP[/3(W// Y.j/^ Pj + ^^^f E 



(A.6) 

yields the same trajectory for the converged states even 
if the intermediate stages and the speed of convergence 
greatly differ. One can also add the increments Ap in 
the reservoir as in the molecular simulation: p^*^™ — 
pcdd ^_ ^j^gjj^ Qjjg computes the new Lagrange mul- 
tiplier with Eq. IA.5I and the new local densities with 
Eq. IA.41 checks the constraint, and iterates using Ap = 

T^[PT - a ^'y;'" - (1 - a)pS'=™]- This does not change 
the isotherm when the evolution is adiabatic. This al- 



gorithm seems simpler and possibly closer to the exper- 
imental protocol. However, it becomes problematic in 
the limit of the canonical ensemble (a=l). Therefore, we 
have preferred the first algorithm described above. 

Our calculations were performed with samples of lin- 
ear sizes varying from L = 25 to L 



convergence was assumed when for the 



max 

{i} 



(«-i) 



10^ 



P: 

6 for L 



(") 



< e and 
25 and e 



Jn-l) _ E.pI 



(") 



10 



Vp 

-8 for L 



100 and 
iteration, 

< e with 

= 50 and 



e 

L = 100. In addition, the steps Ap in p were taken as 
small as 10~^ so that most of the avalanches could be re- 
solved (see Ref. [l| for more details about identification 
of avalanches). Enforcing further the convergence crite- 
rions {e.g. with e = 10~^ for L = 50) does not change 
the path followed by the system: adiabatic regime has 
been reached. However, for the smallest system studied, 



16 

with linear size L = 25, diminishing too much Ap could mctastable state with the required density, 
prevent convergence as the system is not able to find a 



